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CHARACTERIZATION OF BAYES PROCEDURES FOR MULTIPLE 
ENDPOINT PROBLEMS AND INADMISSIBILITY OF THE 
STEP-UP PROCEDURE 1 

By Arthur Cohen and Harold B. Sackrowitz 

Rutgers University 

The problem of multiple endpoint testing for k endpoints is treated 
as a 2 k finite action problem. The loss function chosen is a vector loss 
function consisting of two components. The two components lead to 
a vector risk. One component of the vector risk is the false rejec- 
tion rate (FRR), that is, the expected number of false rejections. 
The other component is the false acceptance rate (FAR), that is, 
the expected number of acceptances for which the corresponding null 
hypothesis is false. This loss function is more stringent than the posi- 
tive linear combination loss function of Lehmann [Ann. Math. Statist. 
28 (1957) 1-25] and Cohen and Sackrowitz [Ann. Statist. (2005) 33 
126-144] in the sense that the class of admissible rules is larger for 
this vector risk formulation than for the linear combination risk func- 
tion. In other words, fewer procedures are inadmissible for the vector 
risk formulation. The statistical model assumed is that the vector of 
variables Z is multivariate normal with mean vector fi and known in- 
traclass covariance matrix E. The endpoint hypotheses are Hi: fj,i = 
vs Ki :Ui > 0, i = 1, . . . , k. A characterization of all symmetric Bayes 
procedures and their limits is obtained. The characterization leads 
to a complete class theorem. The complete class theorem is used to 
provide a useful necessary condition for admissibility of a procedure. 
The main result is that the step-up multiple endpoint procedure is 
shown to be inadmissible. 

1. Introduction. Let Z be a k x 1 random vector which is multivariate 
normal with mean vector \x = (/ii, . . . , /x^)' and known covariance matrix S. 
Assume £ is intraclass, that is, all variances are equal to a 2 and all correla- 
tions are equal to p. Consider the k hypothesis testing problems Hi: m = 
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vs Ki : m > 0, % = 1, . . . , k. This represents the problem of multiple endpoint 
testing. We view this problem as a 2 k finite action problem where we can 
decide whether to reject or accept each Hi individually. The loss function 
chosen is a vector loss consisting of two components. The first component is 
the number of false rejections and the second component is the number of 
false acceptances. The corresponding vector risk has one component related 
to the average power of the procedure while the other component is related 
to the average size of the procedure. This will be made precise in Section 2. 

The vector loss function is more stringent than the linear combination loss 
function used in Lehmann (1957) and Cohen and Sackrowitz (CS) (2005) in 
the sense that the class of admissible procedures for the vector loss function 
contains the class of admissible procedures for the linear combination loss 
function. In other words, any procedure shown to be inadmissible for the 
vector loss is inadmissible for the linear combination loss. 

In this paper we offer a characterization of the class of symmetric (permu- 
tation invariant) Bayes procedures. For the normal model, intraclass is the 
most general case of permutation invariance. The characterization leads to 
a useful complete class theorem. The complete class theorem yields a useful 
necessary condition for admissibility of a procedure. Our most important 
result is that the popular step-up multiple endpoint testing procedure is 
inadmissible. 

Step-up procedures are studied in many places including Hochberg and 
Tamhane (1987), Hochberg (1988) and Shaffer (1995). Step-up procedures 
have been studied in connection with procedures that control the false dis- 
covery rate (FDR). See Benjamini and Hochberg (1995), Benjamini and 
Yekutieli (2001) and Sarkar (2002). Six of the eighteen multiple endpoints 
procedures studied by Dudoit, Shaffer and Boldrick (2003) are step-up pro- 
cedures. 

The inadmissibility result for step-up is of practical importance. Further- 
more, the result is somewhat akin to the Stein-type inadmissibility phe- 
nomenon in the following sense: The step-up procedure leads to admissible 
tests for each component individually when p > 0. See CS (2005). Yet in 
the finite action problem if the loss function is the vector loss function of 
this paper or the sum of losses for the component problems, the step-up 
procedure is inadmissible for k>2. 

We note that the characterization of Bayes procedures for finite action 
problems has only been achieved in the past for the case where the action 
space is a subset of the real line. See, for example, Karlin and Rubin (1956), 
Brown, Cohen and Strawderman (1976) and Van Houwelingen and Verbeek 
(1985). Finite action formulations are realistic, practical and important. 

Preliminaries and notation will be given in Section 2. The characterization 
of symmetric Bayes procedures will be given in Section 3. Section 4 contains 
a description of a complete class, a necessary condition for admissibility, and 
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the result that the step-up procedure is inadmissible. In Section 5, for the 
case k = 2 a procedure that is better than step-up is constructed. 

2. Preliminaries. This 2 k finite action problem has actions a = (ai , 02, . . . , a*,)' 
where ai equals or 1 for i = 1, . . . , k. An action where dj = 1 means that Hi 
is rejected, where if ai = 0, Hi is accepted. Thus, for example, a = (1, . . . , 1)' 
means all H are rejected. It will be convenient to define 

T = {u : u = (tii, . . . , Ufc)', Ui = or 1, all i}. 

Note that T can be used to represent the totality of all actions. However, V 
will serve other purposes as well. 

Decision rules 5(-\z) are probability mass functions on T with the in- 
terpretation that 5(a|z) is the conditional probability of action a given z is 
observed. For each z, a nonrandomized decision rule chooses a single element 
of T with probability 1 and assigns all other actions probability 0. 

Let ipi(z) be the probability of rejecting Hi. A decision procedure <5(a|z) 

determines a set of ip^ (z), i = 1, . . . , k, as follows: 

(2.1) ^(z)=£<K a l z )> 

where Ai = {a € T : a has a 1 in the ith position}. Whereas <5(a|z) determines 
ip(z), the reverse is not true. If tp(z) = . . . ,ipk)' is nonrandomized, it 
uniquely determines some <5(ajz). The <5(a|z) determined is also nonran- 
domized. 

The parameter space is O = {/i : /ij > 0, i = 1, . . . , k}. Partition the param- 
eter space Q, into 2 k sets S1 V5 v£T, where O v = {ju: \i = (pLi, /j,2, ■ ■ ■ , 
/^j > if i>j = 1 and /U, = if = 0, i = 1, . . . , £;}. Also let 0« = {/^ : /x € 
ri,yUj = 0} and let fl^ c be the complement of fiW relative to 17. 

A loss function is a function of the action taken and the true state of 
nature. We will study several different loss functions and their relationships. 
For each individual hypothesis Hi the loss function is zero for a correct 
decision, 1 for rejecting Hi when it is true and 1 for accepting Hi when it is 
false. (Note that the ensuing development would also work if 1 is replaced by 
b, b > 0, when a false acceptance is made.) Such a loss function determines 
a risk 

(2.2) fl(o(^i,A*) = (1 " Vi)EMz) + Vi(l - E^i{z)). 

For the finite action problem a sum of the loss functions of the individual 
problems is 

k k 

(2.3) £(a,/x) = ^Oj(l -Vi) + - a i) v ii A*en v - 

i=l i=l 
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The corresponding risk function is Ya=i > which can be expressed as 

(2.4) E M (V>'(l-v) + (l-V0'v). 

The risk function (2.4) is used by Lehmann (1957) and CS (2005). CS (2005) 
also study a vector loss function consisting of the vector of losses for the 
individual component problems. The corresponding vector risk called VRI 
is 

(2.5) 0R(1)(V>1, A*), • • • , #(fc)C0fc >/•*))• 

Another vector loss consisting of two components (Lq,Li) is where 

k k 

(2.6) L (a,/i) = 5^Oi(l-Ui), Li(a,/i) =53(l-ai)«i, /x G f2 v . 

i=l i=l 

The corresponding risk function can be expressed as (Ro(ip, /i), R\(ip, /z)), 
where 

fc fc 

(2.7) i? (^,/^) = E( 1 -^)^( z )' i?i(^,/x) = ]T^£ M (l-<Mz)). 

i=l i=l 

Suppose to is the number of positive //j. Then according to the definition 
of average power given by Benjamini and Hochberg (1995), and noted by 
Shaffer (1995) and Dudoit, Shaffer and Boldrick (2003), R\jm is 1 minus 
the average power. Ro/(k — m) is the average size. Also one may justifiably 
call Rq or Ro/(k — m) the false detection rate or the false rejection rate 
(FRR) and call R\ or R\jra the false acceptance rate (FAR). We call the 
vector risk (Rq,Ri) in (2.7) VRSP since it is related to average size and 
average power. We note that the class of admissible procedures is largest for 
the VRI risk function in (2.5) and smallest for the risk function in (2.4). Yet 
the class of admissible procedures is certainly larger for VRSP than for the 
risk function in (2.4). Thus any procedure which is inadmissible for (2.7) is 
also inadmissible for (2.4). 

In this paper we focus on VRSP. To deal with VRSP we use a device 
utilized by Cohen and Sackrowitz (1984). That is, we introduce an artificial 
but useful problem. Let 9 be a nuisance parameter which takes on the value 
or 1. Next define the one-dimensional loss function 

(2.8) L*(a,0M))=Ma,/z). 

It now follows from Cohen and Sackrowitz (1984) that the class of admissible 
procedures for the problem using (2.6) as a loss function is the same as the 
problem using (2.8) as a loss function but treating 9 as a parameter which 
can either be or 1. Hence we study the problem using (2.8) as the loss 
function. The corresponding risk function will be denoted as R*(ij), (fj.,9)). 
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Now a decision procedure rj)* is Bayes with respect to (w.r.t.) a prior 
distribution ^((i,9) if 

(2.9) EzR*W*, (/i, 0)) = MEttfty, (», 0)). 
The prior distribution is written as 

J ~Ui(/*)(l-i9), if = 1, 

where (3 is the probability that 9 = and £o(a0 is the conditional distribution 
of (jl given 6 = and where (1 — (3) is the probability that 6 = 1 and £i(ju) 
is the conditional distribution of given 6 = 1. We write the density of z 
given (fj,,6), noting that this density is the same regardless of 9. That is, 
f(z\(i, 0) = f(z\(i, 1) = f(z\(i) where 

(2.10) f(z\(l) = (l/(2^)|S| 1 /2) e -(l/2yE- Ze ,'E-V e -(l/2)M'S-V. 

Note the marginal distribution of z is 

(2.11) f(z)= ( /( Z |/x)[/?^o(/x) + (l-/3)dei (M)]- 
The following theorem describes a Bayes procedure. 

Theorem 2.1. Consider the risk function R*(ip,((i, 9)). The Bayes pro- 
cedure w.r.t. £((/,, 9) is ijj* = (tpl , . . . , ip*)' where 

L 0, otherwise. 

(2.12) 

Proof. The risk function R*(ip, ((i,9)) can be written as 

k k 

R*W, ((i, 9)) = (1 - 9)^(1 - vjE^M*)] +e y £v i E lx [\ - ^i(z)] 
i=i i=i 

(2.13) 



i=l 

To find the Bayes procedure we must minimize the expected risk. Using 

(2.13) we see that this amounts to setting ipi(z) = 1 if the posterior expected 
value 

(2.14) E{l-Q-Vi\z}<0, 
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where now and Vi are regarded as random variables with joint prior dis- 
tribution determined by £(jLt,0). The left-hand side of (2.14) reduces to 

(2.15) P{Vi = 0\z} - P{G = l|z}. 

Now (2.15) is 

/(z|ji)[/3d&(Ai) + (l-/3)dei(M)] 



..aw 
(2.16) 

-(1 -0) j^f(z\^)d^)j/f(z). 
We see that (2.14) and (2.16) lead to (2.12). □ 
The step-up procedure we study is as follows: 

Procedure 2.1. Let Zru < Zr 2 ) <■■■< Zr k ) be the order statistics for 
the set (Z±, . .. ,Z k ) and let Cj be a strictly increasing set of critical values: 

(i) If Zn\ < C\ , accept Hru where Hru is the hypothesis corresponding 
to Z(iy Otherwise reject all Hi. 

(ii) If Hru is accepted, accept Hr 2 \ if ^(2) < @2- Otherwise reject Hr 2 \ , ■ ■ ■ , h 
(hi) In general, at stage j, if Zu\ < Cj, accept Hrjy Otherwise reject 

H {j) ,...,H {k) . 

Call the step-up procedure ■0g U (z). The procedure for k = 2 is shown in 
Figure 1. 



3. Characterization of symmetric Bayes procedures. In order to charac- 
terize symmetric Bayes procedures we first recognize that the problem with 
loss function (2.8) is invariant under the following groups of transforma- 
tions: 

(i) G = {g : gz is a permutation of the coordinates of z; i.e., g is a k x k 
permutation matrix}. 

(ii) G = {g : g(fJ>, 9) is a permutation of the coordinates of n while leaving 
6 as is; Le., g = (g °)}. 

(hi) G = {g: g(a) is a permutation of the coordinates of a, i.e., g = g}. 

Since the problem is invariant under the finite group G, it follows from 
Ferguson [(1967), Theorem 3, page 162] that any symmetric Bayes pro- 
cedure is Bayes w.r.t. an invariant prior distribution. Any invariant prior 
distribution (under g) depends only on the maximal invariant parameter 
, • • • j /•£(*;) j 0) ■ This restriction then implies that for Bayes procedures, all 
prior distributions are symmetric in n for each fixed 0. In particular, the 
conditional distributions £o(a*) an d will be permutation invariant. 
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In order to characterize symmetric Bayes procedures we will be using 
(2.12). We first wish to express the integrand in a simplified fashion. Toward 
this end recall the expression (2.10) for f(z\fj,). Since £ is intraclass, that is, 
S = f j 2 (l-/3)/ + pll / , 1 = (1,...,1)', E- 1 = (o 2 (l-p))- 1 (I-GW) where 
G = p/(l + (k — l)p), we can express the numerator of (2.12) as 

e -(i/ 2 )z' S -z f exp(z / S -i /x _ (i/ 2 ) M 'S- V) 

(3-1) U) 

x[/3^ (/*) + (l-/3)d& (/*)]. 

Noting that z'S"V = (o" 2 (l - P))^ 1 ^'^ ~ Gl'zl'/Lt), letting for fixed z, 

(3.2) d£ fl » = exp{-(<x 2 (l - ^^Gl'zl'/x + (l/2)/i'S- 1 / x}^(/^), 

and without loss of generality taking cr 2 (l — p) = 1, we can rewrite (2.12) as 

, if Q(n«|z) = { / e z '^[/3^(/x) + (l-/3)d?M] 



(3.3) # 



<i-p, 

, 0, otherwise. 

Note that in (3.2) we absorb an expression involving z into the prior. This is 
okay since z is fixed and in the development to follow even when z changes 
l'z will remain constant. 

To characterize symmetric Bayes rules it suffices to consider only sample 
points such that z\ < • • • < Z/.. Now we give: 

Lemma 3.1. Fix z and assume Z\ < • • ■ < Z/.. Then 

(3.4) Q{fl^\z) >Q(ft {m) |z), i = l,...,k-l. 
The inequality is strict unless Z\ = Z{+\ . 

Proof. We need only consider the integral in the numerator of (3.3). 
Writing \fidt&{p) + (1 - (5) d^(/x)] as we note 

(3.5) / e z '»dC(»)= f exp ]T z^ + z i+lfM+1 J d? (/i). 

Make the change of variables fa = fa+i, fa+i = fa in (3.5) to find (3.5) is 
equal to 

(3-6) 

- / eX P E ^j'Mj + z iVi dC(ll)- 

JQ( i + 1 ) \ ■ J ■ ■ . , / 
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Thus from (3.5) and (3.6) we have (3.4). Note that the inequality in the 
proof is strict unless Zi = Zi+i . This completes the proof of the lemma. □ 

Theorem 3.2. Let z be such that z\ < ■ ■ ■ < Let r e {0, 1, . . . , k} be 
the element of the set for which Q(Q( r ) |z) > (1 — f3) > Q(^ r+1 ^ |z), where r = 
means Q(£l^ jz) < (1 — (3) for all i = 1, . . . ,k and r = k means Q(£l^ |z) > 
(1 — (3) for alii = 1, ... ,k. Then the Bayes procedure is ipi(z) = 0, i = 1, . . . , r, 
ipi(z) = 1, i = r + l,...,k. 

Proof. Use Theorem 2.1 and Lemma 3.1. □ 

4. Complete class and inadmissibility of step-up. Symmetric Bayes pro- 
cedures and weak * limits of sequences of symmetric Bayes procedures 
against symmetric prior distributions form a complete class of symmetric 
procedures for this problem. See Weiss [(1961), page 81], where he defines 
a weak * limit as follows: Let tp n be a sequence of procedures. Then Vn 
converges to ip if 

lim R(ilt n ,n) = R(il>,n). 

n — >oo 

Another complete class of procedures for this problem is the set of almost 
everywhere (a.e.) nonrandomized procedures. This follows from a result in 
Matthes and Truax (1967) where it is demonstrated that each admissible 
ipi(z) must be nonrandomized a.e. It follows that the nonrandomized sym- 
metric Bayes procedures and their a.e. limits are a complete class of sym- 
metric procedures for this problem. 

We proceed to give a necessary condition for admissibility based on a com- 
plete class. Toward this end let tj be the following partial sums of (z±, . . . , Zk)- 
That is, let tj = J2i=j ^i, j = 1, - - . , k. Let tk+i = and to = — OO- 
LEMMA 4.1. Let S = {t : t k > tk-i — t k > ■ ■ ■ > t\ — t 2 }. Then for j = 

as a function of tj is strictly decreasing while 
t\, . . . , tj-i, tj+i, . . . , tk are held fixed. 

Proof. Note we may write 

(4 1) Q(rt j) It) = — exp ^ =1 ^J ^ ~ 

/nexp(Eti(i»-ii+i)j"i)^i(M) 

For fixed ti,t2, ■ ■ ■ . . . ,t k the numerator is a strictly decreasing 

function of tj (recall Hj-i > 0) while the denominator, being a Schur con- 
vex function of z (it is convex and permutation invariant), is an increasing 
function of tj, tj E S, while all other partial sums are fixed. See Marshall 
and Olkin (1979) for discussion of Schur convex functions. It follows that 
Q(n&\t) then is a decreasing function of tj . This completes the proof of the 
lemma. □ 
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Lemma 4.2. Let j = 2, . . . , k. Let ij)(t) be a symmetric Bayes procedure. 
Then fort G S, ipj(t) is a nondecreasing function oftj while (t\, . . . ,tj-.i,tj+i, . 
are fixed. 

Proof. Note since ij)(t) is a symmetric Bayes procedure it follows from 
the proof of Lemma 4.1 that iftj(t) is nonrandomized for j = 2, . . . , k. Use 
Lemma 4.1 again to conclude that for t G <S, ipj(t) is a nondecreasing func- 
tion of tj while (ti, . . . , tj-i,tj + i, . . . , tk) are fixed. □ 

Theorem 4.3. Let j = 1, . . . , k — 1. Let t/)(t) be a symmetric procedure 
such that there exists a sample point t* G S for which i/jj(t*) = 0. Then a 
necessary condition for ij)(t) to be admissible is that ipj(t) = for all t G S 
such that tj < tj . 

Proof. Recall that symmetric Bayes and a.e. limits of sequences of 
symmetric Bayes procedures are a complete class of symmetric procedures. 
Now Lemma 4.2 implies that every Bayes procedure has the required prop- 
erty. The required property must also hold for any a.e. limit of a sequence 
of symmetric Bayes procedures. To see this let ij) n (t) = (tpi n (t), . . . ,tpkn(t))' 
be a sequence of symmetric Bayes procedures with i/>(t) its a.e. limit. Since 
ipj n (t) is a nondecreasing function it follows that its a.e. limit is also a non- 
decreasing function. This establishes the theorem. □ 

Corollary 4.4. Let ip(z) be a procedure such that there exists a sample 
point z* = (z*, . . . , zl)' with z* k > z* k _ x > ■ ■ ■ > z\ for which ^(z*) = 0. Then 
a necessary condition for ip(z) to be admissible is that ipk( z ) = for all z 
in the set {z : (z* k + z* k _ x )/2 < z k < z* k , z k _ 2 = 4-2' ■ ■ ■ ,Zi = zf,z k + z k _ x = 

Proof. This follows from Theorem 4.3 since tk = Zk, and fixing tk-i, Zk-2, 
is equivalent to fixing t k _i, . . . , t\. □ 

Theorem 4.5. The step-up procedure given in Procedure 2.1 is inad- 
missible. 

Proof. We show that Procedure 2.1 does not satisfy the necessary con- 
dition for admissibility given in Corollary 4.4. Consider the sample point 
z* = (z*, z%, ■ ■ ■ , z k )' where Zj = Cj — e, j = 1, . . . , k, for some e > to be 
chosen. Note since C\ < ■ ■ ■ < z\ < z^ < ■ ■ ■ < z k and also note that 
■j/? S u(z*) = 0. In particular, the last coordinate of ■j/'su( z *) ^ s zer °- Now 
consider the sample point z = (z* , z%, ■ ■ ■ , z k _ 2 , i&k-iiZk) where z k = z k _\ = 
[{Ck + Cfe-i)/2] — e. Notice that for sufficiently small e, z^-i > C^-i, which 
means that Vsu^) = (0,0, ... ,0, 1, 1)'. In fact there is an open interval of 
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z points on the line Zk + Zk-i = z* k + z* k _ x = t* k _ l ,z k _ 2 , ■ ■ ■ , z\ beginning at 
Zk = *ifc_i/2 and ending before Zk = z k such that ip slJ (z) = (0, 0, . . . , 0, 1, 1)'. 
In particular, the last coordinate of t/> su (z) = 1. This represents a violation 
of the necessary condition for admissibility given in Corollary 4.4. 

The result of Theorem 4.5 is, in a sense, akin to the famous inadmissibility 
result of Stein (1956). Stein considered the model Z ~ N(fi,I) and proved 
that if the loss function is the sum of squared errors, then Z is an inadmissible 
estimator of fi when k > 3. This in spite of the fact that each Z{ is admissible 
for \i{ if the loss function is squared error. In our multiple endpoints testing 
problem ip$y = (^sui> tpSUl, ■■■ ,ips\Jk)' 1S such that ^siJi is an admissible 
test of Hi : m = vs Ki : /ij > when the loss function is (0, 1) and p > 0. 
See CS (2005). Yet ipsv ^ s inadmissible as a finite action procedure when 
the loss is the sum of losses of the component problems [or for the vector 
loss (2.6)]. Here the result is true for k>2. □ 



5. A procedure which beats step-up. In the case of k = 2, the step- 
up procedure is shown in Figure 1(a). It is easily seen that the necessary 
condition of Theorem 4.3 is violated when 2C\ < Z\ + Z2 < C\ + Ci. This 
is the shaded strip in Figure 1(b). By making changes in this strip we show 
how to construct a procedure, ip* , that has a vector risk which is less than 
or equal to the risk of step-up for all fi. 

We begin with any fixed t € (2C\,C\ + C2) and consider Z such that 
Z\ + Z2 = t. Without loss of generality let a 2 = 1. We note that the condi- 
tional distribution of Z\ given Z\ + Z2 = t is iV(| + \{pL\ — ^2), |(1 — p)). 



Also (see Figure 2), for the step-up procedure Vsu( z )) when z\ + Z2 
have 



t, we 



</>su( z ) 



(0,1), ifzi<t-C 2 , 

(0,0), ifi-C 2 <zi<Ci, 

(1,1), ifd<*i<t-Ci, 

(0,0), Xt-C 1 <z 1 <C 2 , 

(1,0), ifC 2 < Zl . 



The procedure ip* is constructed as follows. Consider P f _ ll=f _ l2 (t — C\ < Z\ < 
C 2 \Z 1 + Z 2 = t)- P Ml=A12 (| < Z x < t- d\Zx + Z 2 = t). If we let D(t) be this 



difference in conditional probabilities, then 

2C 2 -t 



(5.1) 



D{t) = § 



Next define C* 



C*(t) by setting 



2$ 



t-2d 
V2(l-P) 



(5.2) 



P[-<Z X <C* 



\Zi + Z 2 



/ 2C*-t 



l-p) 



+ 



5 = W)|. 
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2 su © = (4i) 



T su fe) = (l,l) 



J^iS) = (0,0) 



S su ft)= (!.□) 




(a) 



(b) 



2C,-C 




(c) 

Fig. 1. (a) and (b) show the step-up procedure 0su, (c) shows <j>. 



That is, C* is the solution to 



(5.3) $ 



2C*-t \ 1 



Then for + z 2 = t and D(t)>0, let 



1 / 2C 2 -t 



V72(W) 



(0,1), ifzi<t-C7*, 
</>*(z!,z 2 ) = <( (0,0), ift-C7*<Z!<C* 
(1,0), ifC*<«i. 
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z 2 




c* 

t-C* 



Fig. 2. The step-up procedure with the line zi + zi = t. 

On the other hand, if D(t) < 0, let 

f (0,1), if zi< t-C*, 

t!>*{z U Z2) = \ (1,1), ift-C*<Z!<C% 
1(1,0), ifC*<^. 

The resulting procedure is sketched in Figure 1(c). 

Theorem 5.1. The procedure i/'* is better than tp S u for the vector risk 
VRSP. 

Proof. If we let -0 su (z) denote the step-up procedure, then it will be 
shown that the procedure , 0*(z) above is such that 

(5.4) RoW*,») + A*) < RoWsvV) + bRi(tp sv ,»), 

with strict inequality for some /x, for every b > 0. Note that the procedure 
xj)* does not depend on b. This implies that ijj* beats iJj S y for VRSP. 
Using (2.7) we show (5.4) by showing 

(5.5) £ M {(</> su (z) - "0*(z))'(l - (6 + l)v)|Zi + Z 2 = t} > 
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Table 1 

Evaluation of <j>s\j(z), cj>* (z) and W(z;v) when Z\ + 22 = t 



v 











(0, 0) 


(i, o) 


(0, 1) 


(1, 1) 












1 


-(6+l)v 












(1, 1) 


(-6,1) (1,-6) 


(-6, -b) 














W(z;v) 






^su(z)' 


<P*(z)' 


(4) S u(z)'-0*(z))' 


(4>su 


(z)-tf 


»)'(!- 


(6 + l)v) 


— 00 < Z\ < t — C'2 


(0, i) 


(0, 1) 


(0,0) 














t-C 2 <Z!<Cl 


(0, 0) 


(0,1) 


(0,-1) 


-1 


-1 


b 


b 


Ci<zi<t-C* 


(1, 1) 


(0,1) 


(1,0) 


1 


-6 


1 


-b 


t - G* < zi < C* 


(1, 1) 


(0, 0) 


(1,1) 


2 


1-6 


1-6 


-26 


C* < zi < t - Ci 


(1, 1) 


(1,0) 


(0,1) 


1 


1 


-6 


-6 


t-C 1 <z 1 <C 2 


(0, 0) 


(1,0) 


(-1,0) 


-1 


b 


-1 


6 


C*2 < Z\ < 00 


(1,0) 


(1,0) 


(0,0) 















for all /x G O v , v e T, and t <G {2C±, C\ + C2). We will only study the case of 
D(t) > and \t < C* < t - C x as the other cases are similar. Table 1 outlines 
the possible values that ^* anc ^ (V'su ~~ t l ) *)(l + (&+ l) v ) can take on 
for the possible values of z\, z 2 = t — Z\ and v. Also Figure 2 is helpful. We let 
W(z; v) = (^ su (z) - -ip*(z))'(l - (b+ l)v) and study E„{W{Z; v)\Z l + Z 2 = 
t} as a function of fj, for each v € T. Note that fi is in the parameter space 
only when n € Q v . 

Using the values from Table 1, it is easy to check that the definition 
of C* implies E^ 1= ^ 2 {W(Z; v)|Zi + Z 2 = t} = 0, all v G V. For example, say 
v = (0,l)'. Then, as Zi|Zi + Z 2 =* ~ JV(|, when /xi = /i 2 , 



E fll=m {W(Z;v)\Z 1 + Z 2 =t} 

/ 2d-t \ ; 



+ 



t-2C : 



+ (1-6) 



2C 



V2(W) 
2d-t 



i-2Ci 



VW^pj* 

t-2 

<3? 



V2(W) 
/ 2C 2 -t 



2C* -t 



t-2d 
v/2(l-p) 
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= 



as C* is defined by (5.3). 

When hi ^ \i2 the conditional distribution of Z\\Z\ + Z 2 = t is iV(| + 



We further note that 7/ < when fi € rVi) and r] > when fi € rn o). 
The proof can be completed by studying the pattern of sign changes (see 
Table 1) of W((z±,t — z\); v) as a function of z\. It follows from the varia- 
tion diminishing property [Brown, Johnstone and MacGibbon (1981)] of the 
normal distribution that E^{W{Z; v)|Zi + Z 2 = t} > for all [i G 9, v , v € T. 
This completes the proof. □ 
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we derive our results for the vector loss function of (2.6). The significance 
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